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The usual formulation of thermodynamics is based on the additivity of macroscopic systems. 
However, there are numerous examples of macroscopic systems that are not additive, due to the 
long-range character of the interaction among the constituents. We present here an approach in 
which nonadditive systems can be described within a purely thermodynamics formalism. The basic 
concept is to consider a large ensemble of replicas of the system where the standard formulation of 
thermodynamics can be naturally applied and the properties of a single system can be consequently 
inferred. After presenting the approach, we show its implementation in systems where the interaction 
decays as l/r“ in the interparticle distance r, with a smaller than the embedding dimension d, and 
in the Thirring model for gravitational systems. 

PACS numbers: 05.70.-a, 05.20.-y, 05.20.Gg, 64.60.De 

Additivity plays a central role in the formulation of 
thermodynamics. If a system is divided into different 
parts, each part possessing a certain energy, the system 
is said to be additive if the energy due to the inter¬ 
actions between these parts is negligible in comparison 
with the total energy [l|. Due to additivity, the exten¬ 
sive quantities are linear functions of the system size and 
the thermodynamic potentials present always the same 
concavity. Macroscopic systems with short-range inter¬ 
actions are additive. In contrast, the energy due to in¬ 
teractions between different parts of the system cannot 
be neglected if these interactions are long-ranged, caus¬ 
ing the system to be intrinsically nonadditive. This lack 
of additivity has been identified as the source of the un¬ 
usual thermodynamic properties of systems with long- 
range interactions, typically associated with a curvature 
anomaly of the relevant thermodynamic potential. The 
same happens in small systems with short-range interac¬ 
tions, where the range of the interaction is of the order 
of system size Hi. A peculiar thermodynamic prop¬ 
erty such as negative heat capacity is seen as unusual if 
one takes the thermodynamics of additive systems as a 
paradigm. In practice, properties of this kind are found 
in a large variety of systems in nature, ranging from 
atomic 0! to stellar clusters @-0|, so that they have be¬ 
come just other common properties to be considered. 

Although intense research on systems with long-range 
interactions has been carried out during the last years 
from the statistical mechanics point of view, the thermo¬ 
dynamic framework concerning these systems in connec¬ 
tion with nonadditivity has received much less attention. 

This is mainly due to the fact that statistical mechanics 
has to be necessarily contemplated in order to account 
for the microscopic interactions. Can nonadditivity be 


explicitly identified within the thermodynamic formal¬ 
ism, or one merely has to settle for consider it through 
its implicit contribution to the usual thermodynamic po¬ 
tentials? As we shall see, the thermodynamic formalism 
described here shows the clear role played by nonadditiv¬ 
ity, which can be unambiguously determined and quan¬ 
tified. 

The key idea is to convert the problem of the ther¬ 
modynamics of a nonadditive system to the one of the 
thermodynamics of an additive system and there use the 
standard thermodynamic approach. This can be done by 
considering an ensemble of A" independent, equivalent, 
distinguishable systems. This large ensemble of replicas 
of the system not necessarily has to be interpreted as a 
real physical system; it can be seen as a contrived system 
that helps to infer thermodynamic properties of its sin¬ 
gle constituents. Since the ensemble can be as large as 
needed by taking —> oo, it is in fact an additive sys¬ 

tem and, therefore, the standard equilibrium thermody¬ 
namic approach can be applied. An analogous theoretical 
framework was introduced by Hill for small systems [sj. 
A small system, i.e., a system with a small number of 
particles, is not additive; but additivity is recovered, to¬ 
gether with the usual thermodynamics of macroscopic 
systems, when the number of particles in the system goes 
to infinity, provided the range of the short-range interac¬ 
tion becomes negligible with respect to the system size. 
However, the situation we consider here is clearly differ¬ 
ent: we take from the beginning into account long-range 
interacting systems with a large number of particles. In 
such systems, no matter how large, the size of the system 
and the interaction range are comparable, and therefore 
these systems are always intrinsically nonadditive. 

Let us thus consider a system with energy E , entropy 


2 


S, volume V, and N particles. We now introduce an en¬ 
semble of noninteracting replicas of the systems as a con¬ 
struction from which, as we will see, properties of the sys¬ 
tem itself can be inferred. We stress that, as usual in sta¬ 
tistical ensembles, the replicas do not interact with each 
other. The total energy, entropy, volume, and number 
of particles of an ensemble of such systems are given 
by E t = VE , S t = S , V t = JfV, and N t = JVN, 
respectively. The fundamental thermodynamic relation 
for the ensemble takes the form 

dE t =TdS t -PdV t +fuiN t + gd'jr, (1) 

where T is the temperature, P is the pressure exerted 
on the boundary of the systems, and /i is the chemi¬ 
cal potential of a single system. The last term on the 
r.h.s. of Eq. (JT|| is the central ingredient that this ap¬ 
proach incorporates, which accounts for the energy vari¬ 
ation when the number of members of the ensemble 
varies at constant St, V, and Nt. The function S is called 
the replica energy and quantifies the nonadditivity of 
the single systems; it vanishes for additive systems. To 
see this, consider the following situation. Let us make 
a transformation in which an ensemble of replicas, 
each one with N\ particles, entropy Si and volume V\, 
becomes an ensemble of jV 2 replicas, each one with N 2 
particles, entropy S 2 and volume V 2 , under the assump¬ 
tion that, for a given positive £, we have N 2 = iVi/£, 
S 2 = Si/£, V 2 = Vi/f, but JV, = iJV ,[. Clearly, in 
this case dS t = dVt = dN t = 0, so that, from Eq. Q, 
dE t = $ dJf. But in an additive system the energy is 
a linear homogeneous function of the entropy, volume 
and number of particles, i.e. E 2 = E(S 2 . V 2 , A r 2 ) = 
EiSt/^Vt/^Nt/O = EiSuV^NJ/Z = Ei/e, and 
therefore dE t = 0, requiring $ = 0. Thus, we see that 
additivity implies $ = 0. Hence, S ^ 0 implies nonaddi¬ 
tivity. On the other hand, for a nonadditive system the 
energy is not a linear homogeneous function of the en¬ 
tropy, volume and number of particles, and in general we 
will have $ ^ 0. Thinking for example to the case £ = 2, 
this is a direct consequence of the fact that in a non¬ 
additive system the interaction energy between the two 
halves of a macroscopic system is not negligible. Below 
we will show that § is indeed a property of the system 
under consideration. 

Before proceeding, it is important to stress the fol¬ 
lowing point. We are building a purely thermodynamic 
characterization of nonadditive systems, and we have sin¬ 
gled out one thermodynamic quantity that is peculiar 
for this class of systems. However, we must be aware 
that one of the most striking facts in the statistical me¬ 
chanics study of long-range systems, i.e., ensemble in¬ 
equivalence, should produce a correspondence in a ther¬ 
modynamic treatment, since ensemble inequivalence is 
connected to differences in the macroscopic states acces¬ 
sible to the systems when they are isolated or in contact 
with a thermostat [l|. The difference in the accessible 


macrostates translates into a difference in the equation 
of state between an isolated system and a thermostatted 
one, since, e.g., the temperature-energy relation of an iso¬ 
lated system in the range of convex entropy cannot hold 
for a thermalized system [lj. These caveats do not spoil 
the central role of Eq. <m> in the present treatment; one 
should only consider that all the thermodynamic quanti¬ 
ties in the equations, like e.g. T, S, and $ itself, are those 
corresponding to the actual physical conditions, and can 
be different according to whether the system is isolated 
or thermostatted. 

Equation m can be integrated holding all single sys¬ 
tem properties constant, 


Ed^V = TSdPV - PVd^V + ^NdPV+ SdPV, (2) 


yielding E t = TS t — PV t + nN t +$^Y. Thus, for a single 
system one has E = TS — PV + fxN + S’, together with 
the differential relations 

d S=^dE+^dV-^dN, (3) 

dc? = —S dT + V dP — N dp. (4) 


The first is the usual first law of thermodynamics; the 
second follows by requiring that the differentiation of 
E = TS — PV + fiN + S produces the first equation. 
Moreover, Eq. m shows that the well-known Gibbs- 
Duhem equation for additive systems does not hold here, 
and T, P, and /i may become independent due to the 
extra degree of freedom here represented by $. In the 
context of small systems, this independence between T, 
P, and n has been exploited to consider completely open 
liquid-like clusters in a metastable supersaturated gas 
phase 0- In addition, the deviations of small systems 
thermodynamics with respect to that of macroscopic sys¬ 
tems have been shown in fjJ, [l_2] using the grandcanoni- 
cal ensemble. Furthermore, this approach has also been 
used to study the critical behavior of ferromagnets 0] 
by considering an ensemble of physical subdivisions of a 
macroscopic sample; here we always consider replicas of 
the whole system under consideration. 

In practice, to obtain S’ one computes the entropy S of 
a single system, whence computes the equations of state, 
and thus obtains 


S = E - TS + PV - fiN. (5) 

Furthermore, to compute the entropy in a concrete case 
we will employ the microcanonical partition function. 
Based on the considerations made before, we will then 
describe the thermodynamics of a long-range isolated sys¬ 
tem. 

The microcanonical entropy of a single system can be 
obtained from phase space considerations for the whole 
ensemble. Henceforth, the dimension of the embedding 
space is assumed to be d , and the position and momen¬ 
tum of the particle j in the system k are denoted by 
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Pjk £ and qj k £ R d , respectively. Taking into ac¬ 
count that the jY systems are independent and equiv¬ 
alent, the Hamiltonian H t of the ensemble is given by 
T~L t = i Hfc, where each individual Hamiltonian reads 

Hk = YljLi \Pjk\ 2 /{2rn) + W(q lk , •.. ,gjvfe) with m the 
mass of the particles. Here W is the potential energy of 
a single system that contains the long-range interactions. 
In addition, using units where fcs = 1, the total entropy is 
given by S t = In u> t , where oj t = w t (-®t, V t , N t , ./U) is the 
density of states obtained from the phase space of the en¬ 
semble. This density of states must be computed not only 
considering that the total energy is fixed to E t = ^YE, 
but also that the energy of each single system is fixed 
to E. According to this, the energy constraint in phase 
space can be written as p(E t ) = d(E t /PT — Hk)- 

Since the systems are also considered to be distinguish¬ 
able, the microcanonical density of states of the ensemble 
is therefore given by 


p(Et) 






2dN -p c/K 

1 k = W ? 


( 6 ) 


The mean-field potential energy of a single system can 
be written as W = i / n(x)&(x) d d x, where n(x) is the 
number density and <f>(a;) is the potential characterizing 
the long-range interactions at a point x £ M. d in the one- 
particle configuration space of a single system. Likewise, 
the mean-field entropy of a single system takes the form 

S = — J n(x) In [A ^n(x)\ d d x + ^ ^ N, (8) 

where A t = h/y/2irmT is the thermal wavelength. Since 
long-range interactions act locally as an external field, 
the global chemical potential takes the form 0 

/i = p 0 (x) + $(*), (9) 

which, being constant, prevents net fluxes of particles 
through the system. Here po(x) is the local ideal chemi¬ 
cal potential, which depends on n(x) in such a way that 
from © one obtains n(x) = A^ d exp{— [ < f)(a;) — p]/T}. 
Furthermore, multiplying both sides of © by the number 
density and integrating over the volume one obtains 


where h is a constant, d 2div I4 = d d qjk d d pj k , and 

w = oj(E, V, N ) is the density of states of a single system. 
Since the particles are confined to move within the walls 
of their own system, spatial integrations in ([G]) extend 
over the volumes 14 satisfying J V/ d d qjk = V for all j and 
k. Thus, in view of ©, as required, S\ = ,/F In uj = S , 
which highlights the fact that the information concerning 
nonadditivity is contained in the microcanonical entropy 
S' of a single system. 

The contribution of long-range interactions to ([5]) can 
be further concretized if the thermodynamic quantities 
are separated into a part evaluated without long-range 
interactions and the corresponding excess produced by 
these interactions. Assuming that there are no short- 
range interactions, so that the potential energy is only 
due to long-range interactions, we write the entropy and 
energy as S = + S^ and E = E^ + E^ , respec¬ 

tively, with E = W. Here, quantities labeled with (i) 
correspond to the ideal gas contribution, while (e) indi¬ 
cates the excess produced by the long-range interactions. 
Analogously, we write p = p^+p^ and P = Ph)-|_p( e ). 
Using these expressions, from fl&J) one obtains 

£ = W - TS (e) + P (e) V- p (e) N (7) 


since in the absence of long-range interactions TS h) = 
E (*) + P^V — p^ N. If we include also short-range inter¬ 
actions, and not only ideal contributions, the last state¬ 
ment is still true if the splitting to account for the ex¬ 
cess produced by long-range interactions is performed in 
such a way that Eq. J7]) is satisfied. Expression dTJ) for 
the replica energy can be significantly simplified using 
the mean-field approximation in the large N limit; an 
approximation that can be employed for long-range sys¬ 
tems 


14]. This is discussed next. 


pN = T J d d x n{x) In \X^n{x)\ + 2W. (10) 

Once suitable expressions for the entropy and chemical 
potential have been derived, the relation between their 
excess parts can be explicitly written down. Since S W 
and p W can be obtained from © and ([TU]) . respectively, 
by setting <h(a;) = 0, the excess quantities S^ = S — S® 
and p= p — p ^ follow straightforwardly. As a conse¬ 
quence, one has p^N = — TS^ + 2W, and therefore 

£ = -W + P {e) V. (11) 

Equation m for the replica energy is particularly useful 
since it does not involve entropy and chemical potential. 

As examples of systems where £ can be computed 
exactly using the mean-field description, we have sys¬ 
tems where the interactions decay as 1 jr a in the in¬ 
terparticle distance r, with 0 < a < d. In this case, 
the virial theorem states that dNT + aW = dPV and 
hence P^V = aW/d. Thus, the replica energy be¬ 
comes 0 £ = —(1 — a/d)W, and we see that the sys¬ 
tem becomes additive for the marginal case a = d. In 
the case of infinite-range, nondecaying interactions, i.e., 
a = 0, the system is homogeneous and thus one has 
TS = E W -f pWy _ ^ (OjV. This apparent inconsistency 
is solved by realizing that a = 0 means actually no inter¬ 
action at all, since IT is a constant; therefore the vanish¬ 
ing of the potential at infinite distance requires W = 0, 
and then £ = 0. Also, note that for a > d the mean- 
field approximation fails, and therefore one cannot use 
the last expression to infer that £ ^ 0. 

Another example that we will consider is the Thirring 
model [g] in d = 3. This is a solvable model that incor¬ 
porates some of the remarkable features of systems with 
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FIG. 1. Microcanonical reduced replica energy for the 
Thirring model as a function of the reduced energy for differ¬ 
ent values of 77 . In the inset we show II(A, 77 ) = P^V/^N 2 ), 
which is one of the contributions to the reduced replica energy; 
the other contribution, —W/(uN 2 ), decreases for increasing 
A and rj. 
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FIG. 2. Scaled microcanonical entropy per particle for the 
Thirring model as a function of A for different values of the 
scale factor and the reduced volume 77 . The curves get close to 
each other for values of A for which the replica energy becomes 
relatively small. Even for high energies, nonadditivity makes 
impossible a strict scaling regime. 


long-range interactions. In this model, inside the vol¬ 
ume V of the system there is a core of volume V c where 
all particles in its interior interact uniformly with each 
other. Outside the core, the particles do not interact 
and thus behave as a free gas. For large N, the mean 
field potential can be written as $(x) = —2 i/N c dv c (x), 
where v is a constant, N c is the number of particles 
in V c , and 6y : ( x ) = 1 if x G V c and vanishes other¬ 
wise. Thus, introducing the number of free-gas particles 
N g = N — N c , the potential energy of the system reads 
W = —v(N — ./Vg) 2 . Let us also introduce the reduced 
variables A = E/{yN 2 ) and 77 = \u(V/V c — 1). For fixed 
E , V, and N, and according to the saddle-point method, 
the fraction n g (A 1 rj) = N g /N that defines the equilib¬ 
rium states of the system is obtained by solving [ 6 (] 


3(1 n e ) _ 1 

A + (1 — h g ) 2 



— 77 = 0 . 


( 12 ) 


The microcanonical entropy per particle s(A, 77 ) = 

S(E, V, N)/N takes the form 

s( A , V) = \ ln A + (1 - Tig ) 2 - (1 - Tig) ln(l - h g ) 

—n g In fig + h g 77 — In + 1) + so (13) 
= si (A, 77 ) T s 0 , 

with so = In ( cVN 1 / 2 ), being c a constant. Accordingly, 
the temperature reads as T = ^vN A + (1 — n g ) 2 . 

Since the particles outside the core are free, the pres¬ 
sure P at the boundary of the system is clearly given by 
P(V — V c ) = Nn g T. Thus, from (fill) one has 


g 


^ “ WV 2 = ( X “ + T ["s l 1 + e ’O “ !] . ( 14 ) 


where r = T/ (vN). The function i/j is the reduced mi¬ 
crocanonical replica energy. Using ( 1121 ) to express e~ r> 
in m, it is not difficult to see that for n g > 0 and 
A> 3 (1 — hg), i/j approaches W/(uN 2 ) = — (1 — h g ) 2 , 
which, in turn, decreases in modulus for increasing 77 . 
This behavior can be seen in Fig. [T] for different values of 
the reduced volume 77 . Notice that in the plot if) presents 
a jump at a value of A slightly greater than zero for 
77 = 6 ; this is because the model possesses a first-order 
microcanonical phase transition and the replica energy 
contains this information. 

Furthermore, in order to evince the intrinsic nonaddi¬ 
tivity of the system, it is instructive to study how the 
entropy behaves when the thermodynamic variables are 
scaled. With this purpose, we introduce a scale factor 
£ > 1 and a scale transformation that acts on a quan¬ 
tity Q such that Q' = Since E'/(uN' 2 ) = A/£ 
and V' = £V, the dimensionless parameters A and 77 
transform according to A —> A = A/£ and 77 —>■ 77 = 
ln [£ (e 7 ? + 1) - 1], We can write S'(E', V', N')/N' = 
si (A, ijj+ln (cVAT 1 / 2 ) + | ln£, where the last term comes 
from the scaling of the volume and the number of parti¬ 
cles in so- The fraction fi g now is obtained from (ED with 
A and 77 in the place of A and 77 , respectively. Moreover, 
the entropy per particle of the fully scaled system can be 
expressed as a function of the parameters of the system 
without scaling by defining s(A, 77 ,^) = si(A, 77 ) + | ln£, 
where we have subtracted sq for convenience since it does 
not depend on £. The scaled entropy per particle s is 
shown in Fig. [2j Due to the nonadditivity, it is clearly 
seen that the entropy strongly depends on the scale trans¬ 
formation, as expected. As the energy increases, how¬ 
ever, in this case the system becomes, so to speak, more 
additive since the curves tend to run together, although 
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they never touch each other. If the entropy were a linear 
homogeneous function of E, V, and N, all curves would 
collapse into a single one. The interesting fact is that the 
nonadditivity becomes less noticeable when the replica 
energy is relatively small, as should be expected since 
S = 0 for additive systems. 

While in a statistical mechanics formulation the non¬ 
additivity is naturally codified in the Boltzmann-Gibbs 
microcanonical entropy, as well as in the corresponding 
free energy for other ensembles, we have seen here that, 
in the thermodynamic treatment, nonadditivity emerges 
through an additional degree of freedom, the replica en¬ 
ergy S. According to the differences between isolated 
systems and systems in contact with a thermostat em¬ 
phasized before, we expect that the replica energy will 
depend on the physical situation under consideration. 
Nevertheless, in a long-range systems it will always be 
different from zero. In conclusion, we have shown that 
nonadditive systems can be treated in the standard equi¬ 
librium thermodynamic framework if it is properly for¬ 
mulated. 
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